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ABSTRACT 

We study a combined sample of 264 star-forming, 51 composite, and 73 active 
galaxies using optical spectra from SDSS and mid-infrared (mid-lR) spectra from 
the Spitzer Infrared Spectrograph. We examine optical and mid-IR spectroscopic 
diagnostics that probe the amount of star formation and relative energetic con- 
tributions from star formation and an active galactic nucleus (AGN). Overall we 
find good agreement between optical and mid-IR diagnostics. Misclassihcations 
of galaxies based on the SDSS spectra are rare despite the presence of dust ob- 
scuration. The luminosity of the [Nell] 12.8 pm emission-line is well correlated 
with the star formation rate (SFR) measured from the SDSS spectra, and this 
holds for the star forming, composite, and AGN-dominated systems. AGN show 
a clear excess of [Neill] 15.6 pm emission relative to star forming and composite 
systems. We find good qualitative agreement between various parameters that 
probe the relative contributions of the AGN and star formation, including: the 
mid-IR spectral slope, the ratio of the [NeV] 14.3 pm to [Nell] pm 12.8 fluxes, 
the equivalent widths of the 7.7, 11.3, and 17 pm PAH features, and the optical 
“D” parameter which measures the distance a source lies from the locus of star 
forming galaxies in the optical BPT emission-line diagnostic diagram. We also 
consider the behavior of the three individual PAH features by examining how 
their flux ratios depend upon the degree of AGN-dominance. We find that the 
PAH 11.3 pm feature is significantly suppressed in the most AGN-dominated 
systems. 


S. M. LaMassa is now at Yale University. 
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Introduction 


It has been established that a strong link exists between supermassive black holes 


(SMBHs) and the galaxies in which they live (IFerrarese fe Merrittl l2000l; iGebhardt et al. 


2000l; iTremaine et al.l 120021) and betwee n accreting SMBHs (or active galactic nuclei, AGN) 


and star formation in particular (e.g. iKauffmamn et al.l 120031: ICid Fernandes et al.l l200ll) . 


Clues to this link can be revealed by detailed analysis of the spectra of galaxies and AGN. 
Star formation processes imprint their signatures on the spectrum. Stellar populations of 
various ages produce strong continuum emission in the ultraviolet (UV) through near-IR 
region, while gas photo-ionized by the hot young stars produces prominent emission-lines 
in the UV, optical, and IR. The accretion disk surrounding the supermassive black hole in 
an AGN prod uces its own strong UV and optical continuum emission (the “big blue bump 
(jShieldsI 197811 ). This emission can swamp that of the stellar population of the host galaxy 
when there is a direct view of the accretion disk (i.e. in Type 1 AGN). Dust in the interstellar 
medium of the host galaxy will absorb optical and UV light from both stars and the AGN 
and re -radiate this in t he IR. In AGN, this re processing also occurs in the dusty obscuring 
torus flAntonucci 1993i ; lUrr v _& Padovanilll995l ). In type 2 AGN, the direct view of the AGN 
is blocked by the obscuring torus, and the presence of the AGN can be inferred by prominent 
optical and infrared emission lines from highly ionized gas located beyond the torus and the 
strong mid-IR continuum emission from the torus. In such objects the signatures of host 
galaxys young stars and star formation are also present in optical and mid-IR spectra. These 
type 2 AGN are therefore ideal laboratories for studying the connection between AGN and 
star formation, provided that the tracers of this activity can be disentangled. 

Multi-wavelength star formation indic ators h ave been extensively studied in samples 
of quiescent star forming galaxies (see Kennic uttlll998l . for a review). Star formation rates 
(SFRs) have been calibrated in the optical based on Ha emission, which results from recom- 
bination following ionization due to energetic photons from massive O and B stars, and in the 
UV based on the continuum which provides a window into emission from young stars. Dust 
in the galaxy reprocesses optical and UV photons and re-radiates in the infrared (IR), thus 
revealing obscured star formation. Other IR indicators include emission from fi ne structu re 
lines excited by starburst activity, such as [Nell] 12.81/im and [Neill] 15.56pm (jHo fe Keto 
20071 ). as well as spectral emission features associated with polycyclic aromatic hydrocar- 
bons (PAHs), si nce t hese macro -molecules/small-grains can be ionized by single optical or 
UV photons (ILi fc Drainel 120021) and are prominent in IR spectra of star-forming galaxies 
( Smith et al 2007bl ). 


Such host galaxy star formation diagnostics can be contaminated by the presence of an 
AGN. Photons emitted by the accretion disk enhance the optical and UV emission used to 
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derive SFRs. The obscuring medium enshrouding the central engine contributes significantly 
to the in frared em ission, constituting approximately 20% of the AGN bolometric luminosity 


flSpinoglio fe Malkan 1989t) . Disentangling the contribution from the AGN vs. star formation 
then becomes a necessary intermediate step in using continuum IR emission as a tracer of star 
formation in active galaxies, and such a process can be plagued with uncertainties. Previous 
studies have indicated that emission from PAHs tend to be suppressed in galaxies hosting 


~ — 

the AGN accretion disk destroying PAHs (e.g. 

Voit 

1992|; 

O’Dowd et al. 

2009 

2010; 

LaMassa et al. 

2010; 

Diamond-Stanic & Rieke 

2010 

)• 


A complementary view of the AGN and starburst connection is to study spectroscopic 
signatures that parameterize the relative importance of these two processes. The optical 
BPT diagram (a plot of [OIII]A5007/H/3 vs. [NII]A6583/Ha:) provides a useful diagnostic 
to differentiate between star-forming galaxies, composite systems (galaxies with comparable 
amounts of star- formatio n and AGN activity) and Type 2, or obscured, AGN in the local 
universe ( Baldw in et al. 119811 1. The optical “D” parameter is the distance a source lies 


from the locus of star forming ga laxies o n the_ BPT diagram: a higher D value indicates 


greater AGN dominance flKauffmann et al.l 2003|) . Ratios of IR fine structure lines also 


parameterize the ionizatio n field hardness. For instance, [QIVl26/im and [NeVll4.32qm are 


primarily ionized by A GN flRigbv et al.ll2009l: iDiai non d-S tame et alJl2009[ lArmus et al.ll2 007; 
Goulding et al. 2010) whereas [Nell] 12.81/iin is excited by star formation ( Ho fc Ketol 2007 1. 


The ratio of these lines can then indicate the relative importance of these two processes 

( Genzel et al. 

1998; 

LaMassa et al. 

2010; 

Pereira-Santaella et al. 

2010). The mid-infrared 


(MIR) spectral slope, a 2 o- 30 /im, steepens as the amount of emission from cold dust heated 
by stars increases relative to the amount of emission from hotter dust heated by the AGN. It 
is thus anoth er potentially useful tool to assess the rel ative amount of star formation in active 
systems (e.g. Buchanan et al. 2006 : Deo et al. 2009 1. The equivalent width (EW) of PAH 
grains is another diagnostic to probe the relative amount of star formation to AGN activity. 
An anti-correlation between ionization field hardness and the ratio of PAH features at 7.7/im 
and 1 1.3/iin has been reported and is interpreted as due to an increasing contribution to the 
mid-IR continuum emission by AGN-heated dust. The selective destruction of smaller PAH 


grains may also play a role (e.g. ISmith et all l20Q7bl : lO’Dowd et al.ll2009t iTrever et al.l 12010 


Wu et al. 201 Qj), and so the flux ratios of the different PAH features can then also potentially 
trace the relative contribution of AGN and star formation in active galaxies. 


In this paper, we combine samples of star-forming galaxies, obscured AGN (Seyfert 
2 galaxies, or Sy2s) and composite systems to study the interplay between AGN activity 
and star formation. Using quiescent star forming galaxies as a base-line, we investigate the 
effects of AGN activity on the following star formation diagnostics: optically derived SFRs 
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from the Sloan Digital Sky Survey (IBrinchmann et al.l 120041). the lumi nosities of the IR fine 
structure lines [Nell] 12.81/jm and [Neill] 15.56/jm flHo fc Keto 2007) and the luminosities 
of the p olycyclic aromatic hydrocarbons (PAHs) at 7.7/im , 11.3/rm and 17/un flSmith et al 


2007bl: iFarrah et al.l 120071: iDiamond-Stanic fc Riekel 120101) . Which of these star formation 
proxies agree the best among star-forming galaxies and AGN and are thus least affected by 
the presence of an AGN? 


We expand upon the work of [LaMassa et al . (120101) to test the agreement among di 


agnostics that parameterize the relative contributions of AGN activity and star formation, 


(Smith et al 2007b; O’Dowd et al. 

2009I 

Trever et al. 2010: E 

)iamonc 

-Stanic & Riekel 2010). 

ratios of IR fine structure lines ( Genzel et al. 

1998: Armus et a 

2007; 

Pereira-Santaella et al. 

2010). mid-IR spectral slope (Wn-qn,,™ ! 

3uchanan et al. 20061: 

Deo et al. 2009) and the opt.i- 

cal D parameter (Kauffmann et al. 

2003). In 

LaMassa et al. 

201 o|) we found that Sv2s with 


stronger PAH emission tend to have a softer ionization held and that the MIR spectral slope 
was well correlated with PAH EW. By expanding this parameter space into the regime of 
quiescent star forming galaxies, we test if PAH features and the mid-IR spectral slope are 
dependent on the hardness of the radiation held, regardless of its source, or if ionization by 
an AGN is necessary to appreciably affected observed features. Finally, we use the results of 
this analysis to empirically decompose the mid-IR (MIR) emission into a star-forming and 
AGN component. 


2. Sample Selection 


To assess the efficacy of star formation diagnostics in active galaxies, we have com- 
bined samples of AGN, quiescent star forming gala xies, and c omposite objects. Sources 


Kewlev et al. 

(2001 

) and 

Kauffmann et al. 

(120031): galaxies above the 

Kewlev et al. 

(2001) 


delineation are c 


Kauffmann et al 


Kauffmann et al 


a ssihed as pure Seyfert 2s, sources in between the Kew lev et ajJ (120011 ) and 
(2003) demarcations are classihed as composites and the sources below the 
(120031) boundary are classihed as quiescent star-forming galaxies. 


As discussed in lLaMassa et al .1 (120091) . we selected an optical sample of Sy2s from the 
Sloan Digital Sky Survey (SDSS). These sources were identified as pure AGN on the basis 
of their [OIII]A5007/H/3 and [NII] A6583/Ha hux rati os which place them above the theo- 
retical starburst demarcation from lKewlev et a l. (120010 on the BPT diagram. This sample, 
complete to an [OIII] hux limit of 10~ 14 erg/s/cm -2 (hereafter “[OIII] sample”), comprises 
20 Sy2s. We complement these Sy2s with a mid-infrared (MIR) selected sample, the orig- 
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inal 12/irn sample Spinoglio & Malkanl (1l989h from the IRAS point source catalog, which 
contains 31 galaxies. Using published or SDSS measured values of [OIII]A5007/H/3 and 
[NII]A6583/Ha, we have classified 27 12/irn sources as Sy2s, 2 as composites (F08572+3915 
and NGC 7590) and 2 were dropped due to lack of published [Nil] /Ha values (F04385-0828 
and NGC 3982). As both [0111] and MIR fluxes are indicators of intri nsic AGN power (e.g 
Spinoglio fe Malkanl Il989t iHeckman et al.l 120051; lLaMassa et al.l 1201011 , these Sy2 samples 
should to first order represent unbiased views of the AGN population. The [OIII] sample 
has optical spectroscopy from SDSS whereas the optical data for the majority of the 12/irn 
sample were based on litera ture va l ues. Both samples have low and high resolution Spitzer 
Infrared Spectrograph (IRS .Houck et al. 2004]) data. 


Our comparison samples of predominantly normal galaxies have SDSS spectroscopy, 
Spitzer IRS data and ultraviolet photometry from the Galaxy Evolution Explorer ( Galex ). 
The Spitzer-SDSS- GALEX Spectroscopic Survey (SSGSS, P.I. Schimi novich) contain s 101 


local (0.03 < z < 0 .21) galaxies overlap ping the Lockman Hole region (lO’Dowd et al.l 12009 
Trever et al.l 12010c lO’Dowd et al.l 1201 lh . This sample is augmented by the Spitzer SDSS 


Spectroscopic Survey (S5, P.I. Schiminovich), increasing the sample size of local (in this 
case, 0.05 < z < 0.1) star-forming galaxies to 300 (Bertincourt et ah, in prep). Both the 
SSGSS and S5 samples have IRS Short-Low (SL) and Long-Low (LL) coverage. The S5 
sample also has IRS Short-High (SH) for all sources, but only 1/3 of the SSGSS sources (a 
sub-sample consisting of the brightest galaxies) have high resolution infrared spectra. Based 
on the BPT designations, there are 26 Sy2s, 49 composites and 264 star-forming galaxies in 
the SSGSS and S5 surveys. 

Combining all four samples, we garner 264 star-forming galaxies, 51 composites, and 73 
AGN for this analysis. Throughout this paper, we have adopted a cosmology of H a = 70 
km/s/Mpc, Vt M = 0.27, and = 0.73. 


3. Data Analysis 
3.1. Optical Data 

Optical spectra for the [OIII], SSGSS and S5 samples were drawn from SDSS, using the 
MPA/JHU galaxy value-added catalog ^] This cat alog includes ancillary derived data, such 
as star formation rates (SFRs, Bri nchmann et ah 2004). which we utilize in the subsequent 
analysis. Only four of the 12/irn selected sources have SDSS spectra (3 Sy2s and 1 composite). 


2 http://www. mpa-garching.mpg.de/SDSS/DR7/ 
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Optical data for the remainder of these sources were culled from the literature, as noted in 
LaMassa et ahl (j2010h . Hence, these 24 Sy2s and 1 composite are omitted from the SDSS 


SFR analysis in Section 4. We emphasize here that all SDSS-derived quantities are based 
on the SDSS spectra alone. 


3.2. Infrared Data Reduction 


m 


The details o f the S yitzer data reduction f or the [OIIIl and 12pm sam ples were presented 


LaMassa et all (120101) . the SSGSS sample in lO’ Dowd et all (120091.12011) and the S5 sampl e 


in Bertincourt et ah, in prep., which follows the same procedures as [O’Dowd et a l.l (120091 ). 
In short, the IRS data were run through the standard Spitzer calibration software. Sky 
subtraction was performed for all low-resolution data and for the high resolution sources that 
had dedicated off source observations. All but 68 observations, from the SSGSS sample, have 
high resolution data, leaving 223 star-forming galaxies, 33 composites and 64 Sy2s with high 
resolution IRS coverage. Only eight sources, from the 12p sample, did not have dedicated 
high resolution background observations, but since we utilize these data to measure emission 
line fluxes, the addition of th e backgro un d to t he underlying continuum has negligible effect 
on our derived results (see lLaMassa et al.ll2010h . S pectra were extract ed with SPICE for the 
[OIII], SSGSS and S5 sources and wit h SMA RT ( Higdon et ah 2004 ) for the staring mode 


12pm observations. CUBISM ([Smith et al. 
spectral mapping observations. 


2007al ) was used to extract spectra from 12pm 


3.3. Emission Line Measurements 

Emission line fluxes used in this analysis ([Nell] 12.81pm, [NeV] 14.32pm and [Neill] 
15.56pm) were measured from the high resolution spectra, specifically the Short-High (SH) 
module. The local continuum around the emission feature was fit with a local zero or first 
order polynomial. This local continuum fit was frozen and the emission feature was then fit 
with a Gaussian. The error was estimated by smoothing the local continuum by the line 
width (or the nominal minimum resolvable line width|f| in the case of non-detections) and 
calculating the RMS around this smoothed continuum. In order for a line to be considered 
detected, it had to be significant at the 5cr level or above, have a line width at least on the 
order of the SH resolution (i.e. full width half max > 0.01 pm) and have a line centroid 


3 width = A„ om /R, where A nom is the nominal wavelength of the emission line and R is the resolution of 
the SH module, ~600 
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near the nominal wavelength of the emission line (i.e. < 0.03/un). The [Nell] and [Neill] 
fluxes and upper limits for the SSGSS, S5, [OIII] and 12 /an samples are listed in Tables |T| 
- m Since we consider [NeV] upper limits in aggregate rather than individually (see Section 
7.1), we list only the [NeV] detections in these tables. 


3.4. PAH Measurements 


We measured the fluxes and equivalent widths (EWs) of the polycyclic aromatic hy- 
drocarbons (PAH) at 7.7/an, 11.3/an and 17/un with PAHFIT, a spectral decomposition 
routine that models dust continuum, starlight con tinuum, unr esol ved emission lines, PAH 


emission features and silicate absorption/emission ((Smith et all2007bl h PAHFIT requires a 


single stitched low-resolution spectrum. The Long-Low (LL) spectra was therefore shifted 


j 20% or lower 


see 


O’Dowd et al. 

2011; 

LaMassa et al. 

2010) 


This choice of normaliza- 


tion was made since the SL and SH apertures sizes are well-matched, better enabling us to 
compare quantities derived from the high and low resolution spectra. 

The routine was run with a mixed extinction model on all sources. We note that in 


LaMassa et al.l (20101 ). PAHFIT was run with a screen extinction, which assumes a fore- 
ground absorbing slab rather than dust mixed with gas, starlight and grains. We have 
re-run PAHFIT on the [OIII] and 12/un sources using the default mixed extinction option 
to be consistent with the SSGSS and S5 samples, though differences in results were slight. 
For these two Sy2 sampl es, we also a dded a dust component at 1000K since this may better 
model strong AGN (jPeo et al . 20091 1. though we note the difference in the fluxes and EWs 
between the default temperature array and this modified version is negligible for a majority 
of the sources and a factor of 2 or lower for most of the remaining sources. We consider 
a PAH feature to be detected if it is present at greater than the 5 a level. In the subse- 
quent analysis, we consider PAH upper limits on the flux, but only include EWs for detected 
sources. PAH measurements are listed in Tables [3] - 0 


4. Parameterizing Star Formation: Effects of AGN Activity 


Optical and IR fluxes have been utilized to calibrate SFRs in quiescent galaxies. Below 
we investigate if emission from the AGN contaminates these parameters and if composite 
systems are more similar to star-forming galaxies or Seyfert 2s. Here we con sid er op tical 
SFRs derived from SDSS using the methodology described in Brinchmaiin et al. (120041 ) . We 


also explore IR star formation proxies, namely the fine structure [Nell] 12.81/im and [Neill] 
15.56/im emission lines as well as the luminosity of the polycyclic aromatic hydrocarbon 
(PAH) features at 7.7, 11.3 and 17 /an. 


4.1. SDSS SFR 


fiber • 


Circumnuclear Star Formation 


For star-forming galaxies, SDSS SFR fn, er is calculated from the Ho luminosity, corrected 
for dust extinction, measured through the 3” SDSS fiber. In these systems, Ha emission 
results from the recombination of ionizing photons from massive O stars, making this flux a 
usefu l parameterization of recent star formation. Us ing a Kroupa IMF flKroupa fe Weidner 
20031 ). Ha luminosity is then converted into a SFR (IBrinchmann et al. 2004j). 


Such a straightforward calculation fails for active systems due to emission powered by 
the AGN contaminating the Ha flux. However, a SFR can be estimated in these systems 
using the strength of the Dn(4000) feature, which is a break in the optical spectrum at 4000A 
that indicates the age of the host galaxy stellar population. The SFR fib er for composites and 
AGN in SDSS is calculated starting with the strength of the Dn(4000) break. This is then 
converted into a SFR by adopting the relationship between SFR^r/M* uher and Dn(4000) 


measured for quiescent systems and then using measured value of M* jn, er (IBrinchmann et al. 
2004h . 


4.2. IR SFR diagnostics: Neon Emission Lines and PAH Luminosities 

In the infrared, [Nell] at 12.81/nn has an ionization potential of 21.56 eV and is therefore 
excited by star formation activity. [Neill] at 15.56/nn has a higher ionization potential than 
[Nel l] at 41 eV, but i t beco mes more dominant in lower- metallicity star forming systems (e.g. 
Wu et al. 2006 1. Ho fe Keto ( 2007 1 have demonstrated that the sum of [Nell] and [Neill] 
correlates better with total infrared luminosity (a reliable proxy of SFR) than just [Nell] as 
this sum better accommodates the low metallicity systems. They subsequently calibrated a 
SFR based on [Nell] + [Neill], 

PAH macro -molecules/small- grains are ionized by single ultraviolet or optical pho- 
tons (Li & Drainc 2002l l and are very prominent in the spectra of star forming galaxies 

SFRs have consequently been calibrated on PAH luminosities (e.g. 


( Smith et all 


Farrah et al. 


2007 b) 


20071) . However, the PAH luminosity may be affected by AGN act i vity, and the 


effect may be different for the diferen t PAH features (e.g. lO’Dowd et al.ll2009t iTrever et ah 


20101; iDiamond-Stanic fe Riekel l2010l) . In this analysis, we consider the sum of the PAH 
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luminosities at 7.7, 11.3 and 1 7 jun (L-zp ah 7.7, 11.3, n^m) ■ In Section 5, we examine each PAH 
feature individually as a function of increasing ionization field strength. 


4.3. Comparison of Optical and IR SFR diagnostics 


We compare the optical and IR SFR diagnostics pair-wise in Figures [0- [3] as both scatter 
plots (left hand panels) and histograms of the ratios of the two proxies under consideration 
(right hand panels) to search for systematic offsets among galaxy typeB In these plots, 
blue represents star-forming galaxies, green denotes composites and red indicates AGN. 
To take into account upper lim i ts, we used survi val analysis (ASURV, iFeigelson & Nelson 
1985 ; Isobe and Feigelson 199p|; Isobe et al 1986 1 to determine the best-fit trends to the 


star-forming galaxies and Sy2s (shown as the blue dotted-dashed lines and red dashed lines, 
respectively, in the left hand panels of Figures Q]-|3]), to calculate the mean ratio of SFR prox- 
ies for each galaxy type (Table ED and to quantify the significance of agreement/disagreement 
of the distributions of the SFR proxy ratios among star-forming galaxies, composites and 
AGN (Table ITOT) . We use the Gehans Generalized Wilcoxon Test to calculate the probabil- 
ity that the distributions of SFR proxy ratios among the sub-samples being considered are 
drawn from the same parent population. A probability value (P) < 0.05 indicates that we 
can reject the null hypothesis and that the distributions differ at a statistically significant 
level. We note that not all sources are included in each calculation, either due to absence of 
data (i.e. 25 12/irn sources have no SDSS coverage) or lack of significant detections in the 
diagnostics being considered. The number of sources used for each calculation are listed in 
Tables [9] - [TO] 

As shown in Figure Q] and Table [TUI the [Nell] luminosity in relation to SDSS SFR f iber 
is consistent among star-forming galaxies, composites and AGN. However, the sum of [Nell] 
and [Neill] does show signs of AGN contamination as these sources are systematically offset 
to higher L [Neii] + L [Neiii] values with respect to star-forming galaxies, as seen in Figure 
El Indeed, the distribution of Li NeII ] + L m e nn normalized by SFR f iber is significantly dif- 
ferent between star-forming galaxies and AGN (Table [TUI) . This result is due to the harder 


dicator in active systems. This result is consistent with 

previous works (e.g. 

Deo et al. 

2007; 

Tommasin et al. 

2008; 

Melendez et al. 2008: 

Wu et al. 

2010; 

Pereira-Santaella et al. 

2010) 


4 In these plots and subsequent plots of PAH fluxes and EWs, the error bars are smaller than or on the 
order of the symbol size and therefore not plotted. These relatively small errors are due to the stringent 5tr 
threshold we required for detection. 
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that demonstrated that [NeIII]/[NeII] is correlated with AGN activity and thus a tracer of 
ionization held hardness. Interestingly, the composite galaxies are not affected: their I ^m e ii] 
+ ^ [Neii i] sum differs significantly from AGN, but not from star-forming galaxies. Hence, 
the SFR calibrated on Ljjve//] and L [Nein] is reliable for star-forming and composite galaxies, 
but not AGN. Consequently, to test the efficacy of PAH luminosities in parameterizing star 
formation in AGN, we compare PAH luminosities to SDSS SFR fiber and L[ Ne n\ and not the 
sum of L[jv e jj] and L [Neiny 

Both [Nell] and PAH luminosities are upper limits for at least a handful of sources. 
Since ASURV is unable to accommodate doubly censored data points, we calculated the ratio 
means and performed two sample tests using upper limits on [Nell] and PAH luminosities 
separately. We report the results of both analyses in Tables [9] and [TO] including the number 
of sources for each calculation. We find that the distribution of PAH luminosities relative 
to SFR /ifcer is consistent among star-forming galaxies, composites and AGN, though the 
slope of the relation for Sy2s differs markedly compared with star-forming galaxies (Figure 
[3]). The difference in PAH emission becomes more pronounced when compared to ~L[Neii]i 
regardless of whether L^ PA h 7 . 7 , 11 . 3 , 17 ^ or L [jv e //] is treated as the upper limit (see Table [TUI) . 
As Table [U] demonstrates, there is a decrease of about 32% in the mean in the ratio of PAH 
to [Nell] luminosity in AGN compared to both the quiescent star-forming galaxies and the 
composites. At face value we can not say whether this is caused by supression of the PAH 
emission or enhancement in the [Nell] emission. We will return to this issue later, when we 
explore the effects of ionization held hardness on the relative strengths of PAH features at 
7.7, 11.3 and 17 pm individually. 


5. Relative Importance of AGN to Star Formation Activity 


Expanding on the work of lLaMassa et al.l (120101 ). we now investigate the agreement 
among multi-wavelength diagnostics of ionization held hardness to determine where “nor- 
mal” star-forming galaxies live in this parameter space. As AGN activity weakens and star 
formation becomes more dominant, do the trends found for AGN smoothly extend into the 
regime of quiescent galaxies? 


As discussed in 

LaMassa et al. 

O 

1 — 1 

O 

CN 

), the [OIV]26pm/[Ne 

'll 12. 81 wm hux ratio is a reli- 

able tracer of the incident radiation held hardness ( Genzcl et al. 

1998|; 

Pereira-Santaella et al. 


2010' 

since [OIV] is primarily ionized by AGN activity ( 

Riiibv et al. 20091: 

Diamond-Stanic et al. 

2009 

Melendez et al. 

2008; 

LaMassa et al. 

2010 

) whereas [Nell] 12.81pm is excited by star 

formation ( 

Ho & Ketc 

2007 

). However, we do not have Long-High (LH) coverage for the 


SSGSS and S5 samples, and thus lack high resolution data for the [OIV] line at 26pm. 
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We therefore use [NeVl/ [Neill as a diagnostic of the relative imp ortance of AGN versus 


starburst activitv (Armus et al. 

2007: 

Pereira-Santaella et al. 2010). We compare this ratio 

with other diagnostics studied ii 

1 LaMassa et al. 

(2010). First we use the optical D param- 

eter (D = [log ([N II]/ Ha) + 0.45] 2 + [log([OIII]/ H (3) + 0.5] 2 ), which is the distance a 


source lies on the BPT diagram from the locus of star- forming g alaxies, so that a higher D 
value indicates a greater dominance of AGN activity rfKauffmann et al.l 120031) . Second we 
use O 20 — 30 ,/rn (MIR spectral i ndex) which increases i n the presence o f cold dust associated 
with star formation activity (IBuchanan et al.l [20061: Deo et ah 20091 ). Finally we use the 
PAH equivalent widths (E Ws) w hich are expected to decrease as AGN activity becomes 


more dominant (e.g. iGenzel et al.l Il998l) . We also investigate if the ratios of the fluxes o f 


the different PA H fe atures are affected by AGN activ ity, as found bv lO’Dowd et al.l (120091) . 


Wu et ah (I 2 OIO 1 ) and Diamond-Stanic & Riekg (120101 ). 


5.1. [NeV] 


Due to its hi gh ionization potential (97.1 eV 


SMBH accretion (IGenzel et al.lll998l: lArmus et al 


be an unambiguous signature of AGN activity (IGilli et al. 


, [NeV] has been clai med to cleanly trace 


2007 : Pereira-Santaella et al. 2010 ) and to 


2010 


Gould ing et al. 20 lQl) - Also, 


as it is formed in the narrow line region, it should ideally be unaffected by toroidal obscura- 
tion a nd therefore a tracer of intrinsic AGN flux, s i milar to the [OIIIl 5007A and [OIVl 26/un 


2009 


lines ( Heckman et al.l 12005c iMelendez et al.l 120081; iRigbv et al.l 120091: iDiamond-Stanic et ah 


LaMassa et al.l i2010f) . However, in a majority of our spectra, this feature was not 


significant above the noise. To use [NeV] /[Nell] as a tracer of ionization held hardness, we 
consider the non-detections in aggregate, stacking the spectra among sub-groups of galaxy 
type, rather than employing individual upper limits. For this stacking exercise, we separated 
star-forming galaxies, composites and AGN into two groups each, based on their location in 
the BPT diagram, in order of increasing nuclear activity. Figure [4] illustrates the method we 
used. The l eft hand plot is a BP T diagram for our full sample, with th e unbroken black lines 
showing the iKewlev et al.l (1200 ll) (upper) and iKauffmauu et a, hi (120031) (lower) demarcations 
between star-forming galaxies, composites and AGN. The dashed lines illustrate the addi- 
tional empirical separations we used, approximately tracing the halfway point within each 
region. These lines correspond, in order of lowest to highest, to the following equations: 


Log([OIII\/H(3) = 
Log ([OI 1 1] / H (3) = 


0.61 


Log ([Nil] /Ha) + 0.1 
0.61 

Log([NII]/Ha) - 0.2 


+ 1.3 


+ 1.3 
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Log ([01 II] /HP) 


0.61 

Log ([Nil] /Ha) - 1 


+ 1.19 


We also employed a cut on [Nil] /Ha to exclude low metallicity systems: Log ([Nil] /Ha) 
> -0.6 dex for star-forming galaxies and Log ([NII]/Ha) > -0.2 dex for AGN. The right- 
hand panel of Figur e 0] shows ou r resulting sub-groups d epicted in different colors with the 


Kew lev et al. 020011 ) (upper) and Kauffmann et ah 020031 ) (lower) demarcations overplotted: 


star-forming group 1 - dark blue, star-forming group 2 - light blue, composite group 1 - dark 
green, composite group 2 - light green, AGN group 1 - gold, AGN group 2 - red. 


We stacked together the IRS SH spectra for the sources in each color bin from Figure [4] 
b) not detected in [NeV]. Each spectrum was shifted to the rest-frame, rebinned to a common 
wavelength grid of equal pixel spacing and normalized by the the average MIR flux between 
13.25/mr and 13.75/nn; this window was chosen since it is free of strong emission features. 
Figure 0 shows the stacked SH spectra for each galaxy bin. We used a 3cr clipping on the 
average flux in each pixel to reject outliers that would have otherwise introduced erroneous 
features into our stacked spectra due to noise. [NeV] was detected above the 5cr level in just 
one of these bins, AGN group 2. For the other bins, we obtained [NeV] upper limits using the 
procedure outlined in Section 3. [NeV] /[Nell] was calculated by normalizing the [NeV] upper 
limit/detection by the stacked [Nell] emission feature in each group, which was detected at 
high significance in all cases. In the following analysis, we plot the [NeV] /[Nell] detections 
and stacked [NeV] /[Nell] upper limits against the optical D value, PAH EW (for the sub-set 
of sources where the PAH luminosity was detected above the 5a level) and a 20 - 30 /wn- For the 
stacked [NeV] /[Nell] ratios, we calculate the average D, PAH EW and a20 — 30 pm value in 
each bin. The [NeV] stacked detection is illustrated by the triangle in the subsequent figures. 

Figures [6] - d illustrate the trends between [NeV] /[Nell] and the other diagnostics of 
radiation field hardness. As expected, [NeV] /[Nell] is significantly correlated with D and 
significantly anti-correlated with a 2 o- 3 o^m and PAH EWs, as seen in the correlation coef- 
ficients (p) and associated probabilities that the two quantities are uncorrelated ( P U ncorr ) 
reported in Table [lTJ We note that one galaxy classified as a quiescent star-forming galaxy 
based on its optical (SDSS) spectrum has a [NeV] detection. This suggests that it is har- 
boring a “hidden” AGN, i.e. an AGN embedded in very dusty host galaxy that attenuates 
the optical signatures used to identify AGN. However, we note that this source has other 
IR parameters (i.e. a 2 o- 3 o A im, PAH EWs) more consistent with those of typical star-forming 
galaxies. This may argue against “hidden” AGN interpretation, which may then call into 
question using solely [NeV] detection as an unambiguous signature of AGN activity. 


In any case, we emphasize that the optical and mid-IR classifications of galaxies as 
being dominated by star formation or by a Type 2 AGN agree in the overwhelming majority 
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of cases. This is good news, since it shows that SDSS classifications (which are available 
for roughly one million galaxies) are highly reliable. We do not find a large population of 
optically hidden AGN, as found by Goulding & Alexander (2009). Over half of the [NeV] 
detected sources in their volume limited sample of nearby galaxies are not optically classified 
as AGN. However, they do not include optically classifed LINERS or composites as “optically 
classified AGN,” while we do consider these to contain AGN. The fraction of [NeV] detected 
sources in their sample that are classified optically as quiescent star-forming galaxies is 24% 
(4 of 17), which is still a substantially larger fraction than what we find. 


5.2. 


«20-30/wn 


The m id-IR spectral slope index, « 2 n- 3 n//.m, in creases with the relative amount of star 
formation ([Buchanan et al.l 20061 ; Deo et al. 1200911 . We plot a 2 o- 3 o pm as a function of the 
optical D parameter in Figure [S] which shows that a wide range of a 2 o- 30 /Lcm values are present 
for the star-forming galaxies. Consequently, there is no significant anti-correlation present 
in the full sample, with correlation coefficient (p) of -0.195. However, if we consider just 
the AGN, a significant trend is apparent; p=-0.474. This correspondence likely drives the 
significant anti-correlation found for [NeV] /[Nell] vs a 2 o- 3 o^m- The IR spectral index is thus 
useful in parameterizing the relative contribution of star formation in active galaxies but 
this appears to break down in the quiescent galaxy regime. Consequently, we do not utilize 
« 20 - 30 /im as a tracer of ionization field hardness in the subsequent analysis. 


5.3. PAH EW and Flux Ratios 


Previous studies have indicated that the PAH equivalent widths (EWs) are anti-correlated 
with the presence of increasing AGN activity. This can be due to two effects. First, as the 
AGN becomes increasingly dominant energetically, the emission due to warm dust heated 
by the AGN will rise. Even if the intrinsic PAH emission is unaffected, this increase in the 
mid-IR continuum flux will cause the PAH EWs to decrease. A second possible effect would 
be due to the actu a l dest ruction of dust grains by the hard ionization field from the AGN 
central engine (1Voitlll992h . Our new expanded data set that spans the range from quiescent 
systems to strong AGN provides the opportunity to test these processes over a full range of 
ionization field hardness and to test which features may be most sensitive to AGN activity. 


Figure 0 plots the PAH EW at 7.7 pm, 11.3 pm and 17pm as a function of the optical 
D parameter. The relationship between D and the PAH EWs is significant, with p <-0.65 
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(Table [T2jl . However, if we focus on just the star- forming galaxies and composites, this 
significant anti-correlation disappears: p=-0.227, -0.125 and -0.230 for the 7.7pm, 11.3pm 
and 17 pm features, respectively. Hence, it appears that only in the objects that are optically 
classified as AGN is the PAH EW appreciably affected, suggesting the AGN contributes a 
significant fraction of the mid-IR continuum or perhaps destroys PAHs. 


As noted above, the PAH EW are sensitive to both the relative amount of mid-IR 
continuum from the AGN as well as effects on the PAH emission itself. A complementary 
approach that isolates the effect of the AGN directly on the PAH emission is to consider 
changes in the ratios of the different PAH featur es as the AGN becomes more domin ant. 
Thus, similar to the work of O’Dowd et ah ( 2009 1 and Diamond-Stanic fe Riekel ( 2Qlol) . we 
investigate the relationship between radiation field hardness (parameterized by [NeV]/[NeII] 
and D) and the ratios of the PAH fluxes in Figures [TUI and fill Since the [NeV] line is only 
detected in the optically-classifed AGN (with few exceptions), we compare D to the PAH 
ratios for the AGN subset alone, as well as for the full sample (including star-forming and 
composite systems). Qualitatively, it appears that the distribution of the AGN are different 
in these figures from the star forming and composite objects. For the most AGN-dominated 
sources (D > 1.2 and log(L[jv e y]/L[jv e //]) > —0.5), there is a subset of the population that 
has abnormally low ratios of PAH 11.3 pm to 17 pm and abormally high ratios of PAH 
7.7 /im to 11.3 pm. A similar but weaker effect is seen in the PAH 7.7 /tm to 17 pm 
ratios. These distributions do not appear to take the form of a gradual correlation across 
the whole galaxy population, but rather appear as a threshold effect at high values of D and 
l°g (L[JVeV^]/L[JVe//])- 

Despite this caveat, we have quantified the statistical significances of overall trends that 
may be present. We emphasize that the correlations may be weak, but the trends found can 
be quantified as statistically significant as indicated by the correlation coefficients (p) listed in 
Table fm L P A h7 .7 Mm /Lp ahi7 displays no significant trend with ionization hardness, regard- 
less of the radiation held parameterization used and whether the full sample or AGN only are 
considered. Using [NeV] /[Nell] as a probe of the ionization held, a significant trend becomes 
apparent with L PA Hii. 3 P m/^PAHi 7 u,m (p=-0.564) and ljpAH 7 . 7 P m/^PAHii. 3 P m (p=0.611). How- 
ever, with D as a diagnostic of the incident radiation, the L PA hu .^m/^PAHn^m relationship 
becomes less significant in the full sample (p=-0.364) and even less correlated for the AGN 
sub-sampple (p=-0.247). LpAm^fim/^PAHii^m is not significantly correlated in the full sam- 
ple (p=0.179), but becomes more signhcant in the AGN sample (p=0.431). Taken together, 
these results suggest that the PAH feature at 11.3pm is suppressed in AGN-domin ated sys- 
tems re l ative to the 7.7 and 1 7 pm features. Th is dis agrees with the results of I Smith et al 
(2007b), O’ Dowd et al. ( 2009 1. Wu et al. ( 2010 1 and Diamond-Stanic fe Rieke ( 2Qlol) who 
find that ratio of l^pAH 7 . 7 ^m/^PAHii. 3 fim is smaller in AGN than in star-forming galaxies. 
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The disagreement with the first three studies listed above could result from the inclusion 
of a greater number of high- luminosity AGN-dominated in our analysis. Clearly, AGN affect 
PAH macro-molecu les and this effect co ul d be more pronounced in more luminous AGN. The 
inconsistency with Diamond-Stanic & Rieke 02010!) could be due to their Seyfert sample 
being systematically closer. Perhaps the AGN effect on PAHs are di fferent in the region 
hundr eds of parsecs around the central engine, as probed by the iDiamond-Stanic fe Rieke 
(1201011 analysis, rather than the kiloparsec to galaxy-wide scale that we are investigating 
where there is a greater interplay between host galaxy star formation and AGN activity. 
Different screening effects at different physical scales may play an important role. More 
advanced theoretical work on the effect of the AGN ionization field on PAH macro-molecules 
would help to understand the physical causes of these inconsistencies. 


Because we have combined a number of different samples selected in different ways, our 
AGN systems (mostly drawn from the SDSS [OI1I] and 12/irn samples) are typically at lower 
redshift than most of the star-forming and composite systems (drawn from SSGSS and S5). 
Could any of the results above be influenced by the systematic differences in the physical 
sizes being probed (i.e. smaller scale nuclear emission in closeby AGN versus larger scale 
host galaxy emission for sources further away)? To test this, we impose a redshift cut of 
z>0.03, which corresponds to the lowest redshift source in the SSGSS sample. This cut 
eliminates most of the 12/irn sources and almost half of the [OI1I] sample. However, since 
we do retain a handful of [OIII] -selected AGN which are strong AGN, we are able to test the 
relationship between PAH EWs and ionization field hardness over a wide range of nuclear 
galactic activity. Figure WI\ and I~HT1 illustrate that the trends found for the full sample still 
hold for this pared down sample. The AGN at higher D values are associated with suppressed 
PAH EWs at a statistically significant level even when the host galaxy is more fully sampled 
(see Table IT2l) . The only change in the PAH ratio analysis is the significant anti-correlation 
of LpAHiisixm/^PAHnixm with respect to D in the AGN sub-sample (see Table 1731) . We note 
that though the Lp ahi . 7 iim/^p trend is driven by a handful of AGN, the relationship 
is similar to the parent dataset where the parameter space of high LpAH 7 . 7 pm/^PAHii. 3 pm 
values is more fully sampled. The thresholds in D and [NeV]/[NeII] for changes in the PAH 
ratios noted above are still present. 


Taken together, these results suggest that some care be used when using the PAH 
luminosities to estimate star formation rates in galaxies with dominant AGN. The 11.3 [im 
feature appears to be rather strongly suppressed in the most AGN-dominated systems. 
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6. Disentangling Starburst and AGN Contributions to MIR Emission 


Mid-infrared emission in active galaxies results from both th e toro idal obscuring medium, 
accounting for ~20% of the bolometric AGN flux (ISpinoglio fe Malkan 1989 ). and host galaxy 
star formation. We define L mir as the rest-frame luminosity at 13.5/im averaged over a 3/im 
window: this region is free from strong PAH features and provides a reliable estimate of the 
MIR continuum emission. Using quiescent galaxies as a baseline to estimate the amount of 
MIR emission due to star formation, we can empircally separate the contributions of star- 
burst activity and toroidal emission. To quantify the relative importance of these two effects, 
we solve: 

Lmir = oiL S fr + PLagn 


We set Lrfr = L [Neii] and I jagn = L [oiii],agn, the AGN contribution to the observed 
[OIII] emission. To calculate L[oiii],agn, we use Equation 3 in Wild et ah (120101) . which 
scales the observed [OIII] flux based on the source’s position in the BPT diagram. For 
the AGN where this correction failed (i.e. L[oiii],agn > L [oiii],obs), we set L[oiii],agn = 
L [om],obs- For the composite and star-forming galaxies where this correction failed (i.e. 
L [oiii], agn < 0); we set L [oiii],agn = 0. A vast majority of the star-forming galaxies had 
L [on i], agn set to 0. 


Using the sources that had [Nell] detections, we calculated a and f3 by using a least 
squares trimmed (LTS) regression, a robust regression technique that removes outliers, on 
the above equation. As each execution of the LTS regression results in different values for 
a and (3, we performed this calculation 1000 times to determine the average values and 
spreads on the distributions. The resulting values, using a 3a clipping mean, are a=89±l 
and j3— 111±7. We note that the difference between these values and using a straight average 
are negligible for a and small for /3 (c.f. 109±15). In Figure [Tj we plot L mir vs. oLsfr + 
(3Lagn (he. 89 x L[jv e /j] + 124 x L [oiii\,agn), with the overplotted line indicating where 
the two quantities are equal. No systematic offsets among galaxy sub-type are present. This 
relationship can therefore be valuable for estimating the MIR flux due to star-formation (or 
conversely from the torus) in active galaxies. 


To test this decomposition, we have plotted axL^ NeI i]/LMiR vs. the PAH EW at 7.7/iin 
in Figure [T5l As we have demonstrated above, this PAH macro- molecule/small-grain is the 
least affected by AGN activity (i.e. not suppressed like the 1 1 . 3/irn feature nor mildly boosted 
like the 1 7 /irn feature). Hence, to first order the EW of the 7. 7 /iin PAH feature is tracing the 
“extra” MIR emission from AGN heated dust. If our decomposition is correct, the relative 
amount of MIR emission due to star formation activity, i.e. a x Lnve/n /Lmir, should increase 
monotonically with increasing EW of the PAH 7.7/irn feature. The dotted dashed lines have 
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unity slope and were chosen to intersect the points where the MIR emission is solely due 
to star formation (e.g. log (ax L[ Ne n]/L M iR) = 0 dex) and fiducial PAH EW values that 
approximately enclose the star-forming galaxy locus. Many of the AGN fall within the region 
demarcated by these two lines. The dashed line is the best-fit line to the AGN sources where 
the slope is held constant at unity. This fit neatly bisects the above region. These results 
suggest that our empirical decomposition works reasonably well in describing the relative 
contributions of star formation and AGN activity to the MIR emission. 


7. Conclusions 

We have explored the infrared and optical parameter space where star-forming galaxies, 
composites and AGN live to analyze diagnostics that parameterize host galaxy star for- 
mation. Using star-forming galaxies as a control sample, we have investigated which star 
formation diagnostics are least affected by AGN activity. We have studied diagnostics that 
trace the interplay between AGN and starburst activity and tested whether a smooth transi- 
tion exists over a range of radiation field hardness. Finally, using the results of this analysis, 
we present an empirical decomposition of the MIR flux into a star-forming and an AGN 
component. Our overall result is that the optical and mid-IR diagnostics of star formation 
and of the relative importance of young stars and the AGN generally agree well. Our more 
specific results are summarized as follows: 

SFR diagnostics The SDSS derived SFR jib er and the [Nell] 12.81pm luminosity agree 
well, and are the most reliable SF proxies we have considered for use in star-forming 
galaxies, composites and AGN. The sum of [Nell] and [Neill] 15.56/nn is systemati- 
cally offset to higher values for AGN due to the active nucleus boosting the [Nelli] 
flux. The aggregate PAH luminosity (Lepah 7 . 7 , 11 . 3 , 17 ^) relative to L[jv e //] is sup- 
pressed in AGN compared with both quiescent and composite galaxies. Though the 
(Lsp^H 7 . 7 ,n. 3 ,i 7 A tm)/SFRjj ber . ratio agrees between star-forming galaxies and Sy2s, anal- 
ysis of PAH flux ratios indicate abnormal behavior in some AGN dominated systems, 
such as suppression of the 11.3pm feature, suggesting diminished PAH emission in 
strong AGN rather than [Nell] enhancement. The disparate slopes of (Lj:pAH 7 . 7 ,ii. 3 ,i 7 nm) /^Kfn, er 
between Sy2s and star-forming galaxies also hint that AGN may contaminate the emis- 
sion from PAH macro-molecules. 

Ionization Field Hardness We have used the optical D parameter, which indicates 


the distance a source l ies from the star-fo rming galaxy locus on the BPT diagram 


Kauffmann et al. 

2003 

). and /L[A 7 -„rr] ( Arrrius et al. 

2007; 

Pereira-Santaella et al. 
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20101 ) as probes of the incident radiation field hardness. The infrared diagnostics that 


parameterize the relative contributions of AGN to star-forming activity that we have 
considered (a' 20 - 30 /im; PAH EWs and PAH flux ratios) only show significant trends 
with ionization field hardness in the AGN population. This result suggests that rather 
than parameterizing incident radiation hardness in general, these IR diagnostics are 
only effective descriptors of ionization field hardness in AGN. 


£* 20 — 30 ^™, In AGN, the spectral slope between 20-30 /irn (measured by the spectral 
index has been shown to steepen in the p resence of cold dust due to star 

formation (IBuchanan et al. 20061; iDeo et al. 2 00 9). For the AGN in our sample, 
0 : 20 - 30 iim is significantly anti-correlated with proxies of ionization field hardness, 
albeit with wide scatter, however this result does not hold when expanded into 
the regime of quiescent star forming galaxies. 


PAH EW PAH EWs are significantly anti-correlated with the increasing hardness of 
radiation held, but this is only true when considering the AGN. When considering 
just the star-forming galaxies, no trend exists between incident radiation hardness 
and the streng t h of t he re sulting PAH f eatur e, which agrees with the findings of 


O’Dowd et al. (120091 ) and Trever et al. (120101 ). This result is consistent with the 


hypothesis that the mid-IR continuum has a significant contribution from dust 
heated by the AGN only in objects where the ionizing radiation held is dominated 
by the AGN. For the 1 1. 3/irn PAH feature, our analysis of the PAH hux ratios 
suggests that destruction of the PAH macro-molecule in strong AGN could also 
contribute to suppression of the PAH 1 1 . 3/irn EW. 


PAH flux ratios We And that the luminosity of PAH 11.3 /mi feature is signifcantly 
suppressed relative to the luminosities of the PAH 7.7 /rm and 17 /j,m features 
and relative to the optical and mid-IR derived star formation rates in the most 
AGN-dominated systems. Thus, some care is required in using the PAH 11. 3/irn 
luminosity as a proxy for the SFR in AGN-dominated systems. 


Composite Systems Galaxies that are optically classihed as composite systems are more 
akin to quiescent star forming galaxies in terms of the mid-IR parameters that trace 
star formation: the smaller relative energetic importance of the AGN in composites 
does not seem to affect [Neill] and PAH luminosities. 


[NeV] as unambiguous AGN signature Due to the high ionization p oten tial __ 
(97.1 eV), i t is poten tially an unequivocal signature of AGN activity (IGilli et al. 


2010 


Goulding et ah 2010). One galaxy that is optically classied as a quiescent star-forming 
system in our sample has a [NeV] detection, yet other IR parameters (i.e. PAH EW 
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values) are normal for a pure star forming galaxy. The nature of this galaxy as a strong 
“hidden” AGN is not clear. 


Disentangling MIR Emission Using multiple linear regression, we fit the relation L mir 
= oLsfr + /3Lagn to our full sample to determine a and j3. We set L sfr = L [Neii] and 
^agn=^[oiii],agn where in L [oiii\,agn, we ha ve subtract ed out the estimated starburst 
contribution to the [OIII] flux using Eq. 3 in Wild et ah ( 2010 1. We find a=89±l and 


/3=111±7. This decomposition can be useful in estimating the MIR emission due to 
the circumnuclear AGN obscuration or conversely host galaxy star formation in AGN. 
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Log (SFR fiber M 0 /yr) 



39.5 40.0 40.5 41.0 41.5 42.0 42.5 43.0 
Log ( L [Neii]/SFR fiber erg/ s/M e yr) 


(a) 


(b) 


Fig. 1. — Left: L[jv e /j] vs. SDSS SFRfn, er . Right: Distribution of log(L[7v e j/]/SDSS SFR/^er)- 
The symbol coding is as follows: S5 - crosses, SSGSS - circles, [OIII] sample - squares, 
1 2 /jrn sample - diamonds. Blue data points represent star-forming galaxies, green indicates 
composites and red denotes AGN. The dashed-dotted line indicates the best fit to the star- 
forming galaxy sub-sample and the red dashed line shows the best fit to the AGN population. 
Active galaxies agree with quiescent in terms of the L[jv e /j] distribution relative to SDSS 
SFR /iber . 
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40.0 40.5 41.0 41.5 42.0 42.5 43.0 
Log ((L [Nell] + L [Nelll] )/SFR fiber erg/s/M e yr) 

(b) 

Fig. 2. — Left: L [Neii\ + ^[Nein] vs. SDSS SFR f iber . Right: Distribution of 

log((L[iVe 77 ]+L[jve 7 / 7 ]/SDSS SFR fiber)- We note that the relatively large number of non- 
detections causes the best-fit relation from survival analysis to lie below the majority of 
plotted upper limits. Nevertheless, the AGN clearly have boosted L[ NeII ] + L[ NeIII j values 
compared with star-forming and composite galaxies, which can be seen by comparing the 
best-fit trend to the Sy2s (red dashed line) with the star-forming galaxies (blue dotted- 
dashed line). Compared with Figure 1, a clear offset is discernable between the best-fit 
relations between Sy2s and quiescent galaxies, suggesting AGN enhances [Neill] emission. 
Color, symbol and line coding same as Figure 1. 




Log (SFR fiber M 0 /yr) 
(a) 


25 



Log (SFR fiber M 0 /yr) 
(a) 



41.0 41.5 42.0 42.5 43.0 43.5 44.0 44.5 
Log (Ljpah 7.7,1 1.3,17/im/SFRfiber er 9/ S /^0 Y r ) 


(b) 




Fig. 3. — Top Left: l^spAH 7 . 7 ,n. 3 ,i 7 nm vs SDSS SFR fn, er . Top Right: Distribution of 
l°g(LsPAi77.7,11.3,17/im/SDSS SFR fiber)- Bottom Left: L S pAP7.7, 1 1.3, 1 7/im vs L [Nell]- Bottom 
Right: Distribution of log(L S pAP 7 . 7 ,n. 3 ,i 7 /im/L[Ar e //]). The sum of the PAH luminosities rela- 
tive the SFR/jfe er is consistent among star-forming galaxies, composites and AGN. However, 
when compared with L[j v e pr], the PAH emission in AGN is suppressed by about 32% on- 
average in relation to the observed PAH emission in star-forming galaxies and composites: 
the mean log(L S pAP 7 . 7 ,n. 3 ,i 7 /im)/L[Ar e //] is 1.50 dex for star-forming galaxies, yet 1.33 dex 
for Sy2s. The best-fit trend to the Sy2s (red dashed line) also illustrates the different re- 
lationship between the PAH luminosities in AGN compared with quiescent galaxies (blue 
dotted-dashed line). Color, symbol and line coding same as Figure 1. 
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Fig . 4. — Left: BPT dia gram for all s ources in our sample. Th e unbroken black lines indicate 
the iKewlev et al.l (120011 ) (upper) and iKauffmann et al.l fl200.l l (lower) demarcations between 
star-forming galaxies, composites and AGN. The dashed black lines illustrate the dividing 
lines we used to separate the sample to stack the short-high IRS spectra. Symbol and color 
coding is same as Figure 1. Right: After removing the low-metallicity systems and using the 
demarcations shown in the left-hand plot, the groups into which our sample were separated 
our indicated by the colors in this plot: star-forming group 1 - dark blue, star-forming group 
2 - light blue, composite group 1 - dark green, composite group 2 - light green, AGN group 
1 - gold, AGN group 2 - red. 


Flux density (10 "erg/s/cm 2 /fj,m) Flux density (10 "erg/s/cm Flux density (10 "erg/s/cm 2 /fi m) 
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SF group 1 SF group 2 






Fig. 5. — Stacked SH spectra around the [NeV] emission line at 14.32 pm for each group of 
galaxies illustrated in Figure 0]b). The dashed line indicates the nominal [NeV] wavelength. 
[NeV] was detected above the 5<r level in AGN group 2 only. 
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Fig. 6. — Left: lj[NeV]/^[Neii] vs. the optical D parameter. Right: L[ NeV ]/l^[Neii] vs. 
« 20 - 30 iim- These quantities are statistically correlated and anti-correlated: p=0.562 and 
p=-0.568, respectively. Color and symbol coding same as Figure [JJ The red triangle indi- 
cates the [NeV] detection for the stacked spectra in AGN group 2. For these and subsequent 
plots with IR emission line and PAH fluxes, error bars are smaller or on the order of the 
symbol size and therefore not plotted. 
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Fig. 7. — a) \j[NeV]/^-‘[Neii] vs. PAH EW at 7.7 pm. b) L[ 7 v e v]/L[jv e /j] vs. PAH EW at 
11.3 prn. c) L[ Ne v]/^[Neii] vs. PAH EW at 17 /irn. These quantities are statistically anti- 
correlated: p=-0.702, p=-0.779 and p=-0.628, respectively. The red triangle indicates the 
[NeV] detection for the stacked spectra in AGN group 2. Color and symbol coding same as 
Figure 1. 
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Fig. 8. — a^o— 30 /wn vs. the optical D parameter. These parameters are correlated for AGN, 
p=-0.474, but not for the full sample, p=-0.195. Color and symbol coding same as Figure [4] 
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Fig. 9. — The optical D parameter vs a) PAH EW at 7.7 /mi, b) PAH EW at 11.3 /mi and 
c) PAH EW at 17 /mi. These quantities are statistically anti-correlated: p=-0.679, p=-0.747 
and p=-0.658, respectively. Color and symbol coding same as Figure 1. 
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Fig. 10.— L[jVey]/L[jVei 7 ] vs. a) L PAH7 , 7lMrn /L PAH17lirn (no significant trend, p=-0.241), b) 
LpAHii. 3 /im/LpAHi 7 /im (significant anti-correlation, p=-0.564) and c) ^PAm.i^m/^PAHii.^m 
(significant correlation, p=0.611). The red triangle indicates the [NeV] detection for the 
stacked spectra in AGN group 2. Color and symbol coding same as Figure 1. 
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Fig. 11 — Optical D parameter vs. a) L PAif7 . 7MTO /L PAm7Mm , b) L PAim . 3Mm /L PAm7MTO , P 
0.364) and c) Am n ^ml^p AHw.^m- Color and symbol coding same as Figure 1. 
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(a) (b) 



(c) 

Fig. 12. — The optical D parameter vs. PAH EWs. Imposing a redshift cut, to test for 
possible systematic differences due to physical sizes probed between closeby and far away 
sources, does not alter the trends seen in the full sample. Color and symbol coding same as 
Figure 1. 
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Fig. 13. — For sources at z> 0.03, optical D parameter vs. PAH ratios shown in Figure fl3l 
Imposing a redshift. cut does not alter the trends seen in the full sample. Color and symbol 
coding same as Figure 1. 


MIR 



a L[Nell] + ^L[ 0 |||] iAGN 


Fig. 14. — L mir vs. oLsfr + /3L agn, where L S fr = L[jv e //] and L^gat = L[ojjj], j 4Gjv- « and 
/3 were found by using a least trimmed squares regression: a=89±l, f3= 111±7. The dashed 
line indicates where the two quantities are equal. Color coding same as Figure 1. 


37 



1 .5 


| . , « ' | ' * ' 1 | * * 

i 1 I 1 r- 

- 


1 .0 


o 



CE 






_l 

0.5 


o 



" s '\ 






r TT l 



! 

X 

X 

X 


- 

Z 

_l 

0.0 


0 

^x x 

” 

CD 

- 0.5 


□ o ' '' <P 


- 

oo 





- 

cn 

O 

- 1.0 


Z *■' o oo 

o, - ' o 



_i 

- 1.5 

!" ^ " 

' ^ o 




- 2.0 

- 



- 


-1 

.0 

- 0.5 0.0 0.5 

1 .0 

1 


Log (PAH EW 7.7yu,m) 

Fig. 15. — Log (axL[ Ne n]/LMiR) vs. Log (PAH EW) at 7.7/irn. The two dotted-dashed lines 
have unity slope and intercept points where the MIR emission is due purely to star formation 
according to our decomposition (i.e. ax L [Neii] = L mir) at two fiducial values of PAH 
7.7/im EW (1.0 dex and 0.5 dex) chosen to bracket the star-forming galaxies. The dashed 
line indicates the best-fit through the AGN, with the slope held constant at unity. The AGN 
generally fall within the region indicated by the lines, which is consistent with a monotonic 
decrease in PAH EW, denoting the “extra” MIR emission due to AGN, correlating with 
the empirial MIR AGN parameterization (i.e. decreasing axL^ e jj]/LMm values indicate 
greater AGN contribution to the MIR luminosity). Color and symbol coding same as Figure 
1 . 


Table 1. Emission Line Luminosities for SSGSS Sample 


Source Number 

RA 

Dec 

z 

Type 

[Nell] 12.81/im 
10 41 erg/s 

[NeV] 14.32/xm 
10 41 erg/s 

[Neill] 15.56/zm 
10 41 erg/s 

33 

160.291 

56.932 

0.185 

AGN 

8.12± 0.70 


3.34± 0.57 

61 

160.952 

58.197 

0.073 

AGN 

3.19± 0.08 


0.30± 0.05 

11 

161.720 

56.252 

0.047 

Comp 

0.17± 0.01 


< 0.06 

12 

162.268 

56.224 

0.072 

Comp 

3.04± 0.09 


0.31± 0.04 

15 

161.758 

56.307 

0.153 

Comp 

4.40± 0.30 


< 1.64 


1 This table is available in its entirety in machine readable format on-line. We list several row here as a 
guide. 


Table 2. Emission Line Luminosities for S5 Sample 


RA 

Dec 

z 

Type 

[Nell] 12.81/mn 
10 41 erg/s 

[NeV] 14.32/jin 
10 41 erg/s 

[Neill] 15.56/jm 
10 41 erg/s 

36.380 

-8.419 

0.055 

AGN 

1.16± 0.05 


1.02± 0.07 

36.435 

-7.868 

0.077 

AGN 

2.72± 0.12 


1.03± 0.08 

132.186 

53.746 

0.083 

AGN 

4.72± 0.14 

1.18± 0.09 

2.54± 0.08 

132.714 

54.823 

0.081 

AGN 

3.17± 0.11 


2.25± 0.07 

132.815 

54.986 

0.083 

AGN 

3.64± 0.16 


1.68± 0.09 


1 This table is available in its entirety in machine readable format on-line. We list several 
row here as a guide. 


Tabic 3. Emission Line Luminosities for [OIII] Sample 


Galaxy 

RA 

Dec 

z 

[Nell] 12.81/im 
10 41 erg/s 

[NeV] 14.32/mr 
10 41 erg/s 

[Neill] 15.56yum 
10 41 erg/s 

NGC 291 

13.375 

-8.768 

0.019 

1.61± 0.04 

0.60± 0.02 

1.12± 0.02 

Mrk 609 

51.356 

-6.144 

0.034 

5.51± 0.12 


1.80± 0.06 

IC 0486 

120.087 

26.614 

0.027 

1.10± 0.10 

0.55± 0.06 

1.04± 0.05 

2MASX J08035923+2345201 

120.997 

23.756 

0.029 


< 0.14 

< 0.42 

2MASX J08244333+2959238 

126.180 

29.990 

0.025 

0.54± 0.06 

0.55± 0.07 

0.70± 0.06 


1 This table is available in its entirety in machine readable format on-line. We list several row here as a guide. 


Table 4. Emission Line Luminosities for 1 2 /mi Sample 


Galaxy 

RA 

Dec 

z 

Type 

[Nell] 12.81/jm 
10 40 erg/s 

[NeV] 14.32/jm 
10 40 erg/s 

[Nelli] 15.56/iin 
10 40 erg/s 

IRAS 00198-7926 

5.473 

-79.169 

0.073 

AGN 

35.65± 1.45 

49.05± 1.53 

60.06± 1.31 

NGC 0424 

17.865 

-38.083 

0.012 

AGN 

3.06± 0.19 

4.64± 0.27 

5.67± 0.17 

NGC 1068 

40.670 

-0.013 

0.004 

AGN 

5.34± 0.28 

10.11± 0.48 

15.02± 0.32 

NGC 1144 

43.790 

-0.178 

0.028 

AGN 

10.79± 0.26 

0.61± 0.10 

2.88± 0.05 

NGC 1320 

51.203 

-3.042 

0.009 

AGN 

0.47± 0.02 

0.59± 0.02 

0.78± 0.02 


1 This table is available in its entirety in machine readable format on-line. We list several row here as a 
guide. 


Table 5. PAH Luminosities and EWs for SSGSS Sample 


Src Num 

Type 

L7. 7/xm 
10 41 erg/s 

EW 7.7/iin 
// m 

Lll.3 fim 

10 41 erg/s 

EW 11.3/mn 
/irn 

L 17 fim 

10 41 erg/s 

EW 17/irn 
/irn 

13 

AGN 

16.59± 0.82 

1.23 

7.61± 0.32 

0.39 

6.12± 1.09 

0.27 

21 

AGN 

30.66± 1.73 

2.95 

10.80± 0.60 

1.43 

6.44± 0.49 

1.09 

23 

AGN 







33 

AGN 

157.18T 6.51 

3.14 

53.40± 1.08 

1.01 

37.98± 3.62 

0.66 

50 

AGN 








1 This table is available in its entirety in machine readable format on-line. We list several row here as a 
guide. 


Table 6 . PAH Luminosities and EWs for S5 Sample 


RA 

Dec 

Type 

L 7.7 [im 
10 41 erg/s 

EW 7.7/un 
/irn 

Lll.3 /im 

10 41 erg/s 

EW 11.3/im 
H m 

Ll7 fim 

10 41 erg/s 

EW 17/im 
/irn 

36.380 

-8.419 

AGN 

24.01± 2.34 

5.74 

7.92± 0.71 

3.62 

15.13± 2.76 

2.28 

36.435 

-7.868 

AGN 

59.59± 4.64 

1.44 

16.41± 0.62 

0.53 

7.79± 1.54 

0.28 

132.186 

53.746 

AGN 

144. 14± 5.69 

4.81 

41.27± 0.81 

1.81 

32.40± 1.82 

1.38 

132.714 

54.823 

AGN 

54.59± 5.05 

4.51 

21.61± 0.71 

1.89 

16.96± 1.69 

1.04 

132.815 

54.986 

AGN 

107. 85± 3.30 

5.05 

29.54± 0.49 

1.62 

19.68± 1.23 

1.17 


1 This table is available in its entirety in machine readable format on-line. We list several row here as a guide. 


Tabic 7. PAH Luminosities and EWs for [OIII] Sample 


Galaxy 

L 7.7 (an 
10 41 erg/s 

EW 7.7/iin 

/irn 

Lll.3 [im 

10 41 erg/s 

EW 11.3/rm 
//rn 

L] 

10 41 erg/s 

EW 17/irn 
//rn 

NGC 0291 

34.69± 0.95 

10.07 

7.44± 0.15 

0.75 

9.20± 0.15 

0.50 

Mrk 0609 



52.35± 0.78 

1.08 

54.58± 0.54 

0.94 

IC 0486 

22.69± 0.47 

0.75 

4.15± 0.19 

0.12 

5.77± 0.12 

0.16 

2MASX J08035923+2345201 



1.49± 0.05 

0.30 

2.10± 0.13 

0.40 

2MASX J08244333+2959238 



2.80± 0.16 

0.04 

5.67± 0.29 

0.10 


1 This table is available in its entirety in machine readable format on-line. We list several row here as a guide. 


Table 8. PAH Luminosities and EWs for 12 /an Sample 


Galaxy 

Type 

L7.7/im 

10 41 erg/s 

EW 7.7/irn 
/irn 

Lll.3/im 

10 41 erg/s 

EW 11.3/an 
/irn 

Ll7 fim 

10 41 erg/s 

EW 17/irn 
/irn 

NGC 0424 

AGN 

10.25± 0.65 

0.09 

4.46± 0.35 

0.06 

8.05± 0.61 

0.19 

NGC 1144 

AGN 

258. 85± 1.37 

2.21 

92.63± 0.99 

1.83 

51.43± 0.73 

1.56 

NGC 1320 

AGN 

9.69± 0.38 

0.52 

2.45± 0.08 

0.16 

1.31± 0.12 

0.11 

NGC 1386 

AGN 

0.84± 0.02 

0.71 

0.16± 0.00 

0.14 

0.26± 0.01 

0.27 

NGC 1667 

AGN 

61.97± 1.14 

3.85 

18.76± 0.22 

2.18 

9.18± 0.13 

1.63 


1 This table is available in its entirety in machine readable format on-line. We list several row here as a 
guide. 


Table 9. Ratios 1 of SFR proxies 


Ratio 

SF 1 2 

N 2 Composites 

N 

AGN 

N 

Comments 

L [Nell] / S F R fiber 

40.971°;! 

223 s 

40.911°;! 

32 

40.921°;! 

39 s 


( L [Nell] + L [N el II ] ) / S F R f ib er 

40.84l°; 4 I 

223 s 

40.7VS;“ 

32 

41 ii+ 0 - 21 

39 


LeP AH 7. 7, 11. 3,17 fim/ SFRy^ eT . 

42.50;<>;“ 

261 3 

uz. jy_ 0 30 

46 

42 49+ 0 - 16 
uz.n a_ 0 36 

35 s 


LEPAff7.7,11.3,17/xm/L[JVe/J] 

1.50 4 

20 7 3 

i cq+o.05 

l.Oy_ 0 .22 

31 

-| qq+0.15 
1 *°°-0.59 

21 

PAH ULs 

L[Ar e //]/LsPAP7.7,11.3,17Atm 

-1.61 4 

141 3 

-1 62 +0 ' 04 
i - uz -0.13 

25 

-1 c;n+ 0 - 02 
i.OU_ 0 .21 

17 

[Nell] ULs 


1 Ratios in log space. Ratios were calculated using survival analysis. Error bars represent 
the deviation from the values at the 75th and 25th percentiles from the mean. 

2 “SF” refers to star-forming galaxies. 

3 Number of sources in each galaxy sub-group used in calculation of mean values. 

4 Upper limit data point changed to a detection for Kaplan-Meier computation, which biases 
the mean estimate. 

4 Error in percentile calculation, so error bars are not reported. 


Table 10. Two-Sample Tests on SFR Proxies between Galaxy Sub-Samples 1 


Ratio 

SF vs. Composites 
P 

N 2 

Composites vs. AGN 
P 

N 

SF vs. AGN 
P 

N 

Comments 3 

L[7Ve//] / SFRf iber 

0.127 

255 

0.598 

71 

0.624 

223 


(L[iVe//] +L [Nell I] ) / SFR f iber 

0.175 

255 

0.0005 

71 

<lxl0 -4 

262 


LsPAi/7.7,11.3,17//ra/SFRj^ er . 

0.331 

307 

0.365 

81 

0.758 

296 


LuPAiT7.7,11.3,17/im/L[jVeII] 

0.213 

238 

0.0002 

77 

0.003 

253 

PAH ULs 

L[JVe//]/LspA/f7.7,11.3,17Atm 

0.945 

166 

0.006 

42 

0.0004 

158 

[Nell] ULs 


la SF” refers to star- forming galaxies. 

2 Number of sources in two-sample tests. 

3 We report the survival analysis probability from Gehan’s Generalized Wilcoxon Test, permutation variance. For 
the PAff and [Nell] ratios, survival analysis was performed twice: once to account for [Nell] upper limits and once to 
accommodate the PAH luminosity upper limits. A P value <0.05 indicates that the two sub-populations differ at a 
statistically significant level. 
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Table 11. Correlation between Diagnostics of Relative Importance of AGN Activity to 

Star formation 1 


Diagnostic 

P 2 

p 3 

- 1 uncorr 

L [NeV] /L [Nel /] VS. D 

0.562 

4xl0~ 4 

^[NeV]/^[NeII] VS. a20-30/im 

-0.568 

6xl0 -4 

L[NeV]/L[Neii] vs. PAH EW 7.7 pm 

-0.702 

lxl0“ 4 

L [ NeV] / L [ Nell] vs. PAH EW 11.3 /till 

-0.779 

< 1 x 10~ 4 

L[NeV]/L[NeII] vs. PAH EW 17 /JUl 

-0.628 

2x 10~ 4 


1 Used ASURV to account for upper limits in L[jv e y]- 

2 p is the correlation coefficient from which V uncorr is de- 

rived. 

Probability the two quantities are not correlated. A 
Puncorr value <0.05 indicates a significant correlation. 
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Table 12. Relationship between PAH EWs and Ionization Field Hardness 


Diagnostic 

P 1 

N 1 2 

D vs. PAH EW 7.7 pm 

-0.679 

344 

D vs. PAH EW 11.3 pm 

-0.747 

368 

D vs. PAH EW 17 /irn 

-0.658 

281 

Star-forming Galaxies and Composites Sub-sample 

D vs. PAH EW 7.7 pm 

-0.227 

291 

D vs. PAH EW 11.3 pm 

-0.125 

303 

D vs. PAH EW 17 pm 

-0.230 

221 

Sources with z > 0.03 

D vs. PAH EW 7.7 pm 

-0.566 

318 

D vs. PAH EW 11.3 pm 

-0.653 

336 

D vs. PAH EW 17 pm 

-0.549 

249 


1 p is the correlation coefficient. Higher \p\ values 
indicate more significant (anti-) correlations. 


2 Number of sources. 
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Table 13. Relationship between PAH Ratios and Ionization Field Hardness 


Diagnostic 

P 1 

N 2 

[NeV]/[NeII] 3 vs.: 

I J PAP7.7/im/LpAP17/im 

-0.241 

33 

LpAP11.3/im/LpAP17/im 

-0.564 

37 

L p AP7. 7/jm /■ L P API 1 . 3/im 

0.611 

34 

D 4 vs.: 

LpAP7.7/im/LpAP17/im 

-0.193 

266 

L p AP11 . 3/im / LpAPl 7 Mm 

-0.364 

278 

LpAP7.7/im/LpAP11.3/im 

0.179 

340 


AGN sub-sample, D 4 vs: 


L p AP7. 7/im / L p API 7/im 

0.082 

47 

LpAP11.3/im/LpAP17/im 

-0.247 

58 

LpAP7.7/tm/LpAP11.3/im 

0.431 

51 

Sources with z > 0.03, 

D 4 vs.: 


LpAP7.7/jm/LpAP17/jm 

-0.283 

240 

LpAP11.3/im/LpAP17/im 

-0.356 

246 

L PAP7. 7/im / L p API 1 . 3/im 

0.084 

314 

AGN with z> 0.03, D 4 

vs: 


LpAP7.7/tm/LpAP17/im 

-0.208 

22 

L p AH 1 1 . 3/im / LpAP 1 7/i m 

-0.554 

27 

LpAP7.7/im/LpAP11.3/im 

0.590 

26 
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l p is the correlation coefficient. 
Higher \p\ values indicate more signifi- 
cant (anti-) correlations. 

2 Number of sources. 

3 p calculated from survival analysis. 

4 p calculated from linear regression. 


